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Abstract 

The tunneling rates for scalar fields with quartic potentials in de Sitter 
space in the limit of no gravitational back reaction are calculated numeri- 
cally and the results are fitted by analytic formulae. 



1 Introduction 



The decay of a false vacuum and/or creation of topological defects through quan- 
tum tunneling in the early Universe may play an important role in determining 
the subsequent structure of the Universe. These phenomena are particularly im- 
portant in the context of inflationary cosmology. The nucleation of true vacuum 
bubbles associated with false vacuum decay is incorporated in several viable mod- 
els of inflation as a mechanism to terminate the inflationary expansion of the Uni- 
verse 0. Quantities such as the efficiency of thermalization by bubble collisions 
or the degree of spatial inhomegeneities crucially depend on the nucleation rate. 
Also, topological defects such as domain walls, strings, or monopoles may be cre- 
ated during the inflationary stage of the Universe . In order to quantify their 
effects on the structure of the Universe, it is necessary to know their formation 
rates. 

The study of false vacuum decay was initiated by Voloshin, Kobzarev, and 
Okun ||, and subsequently developed by Coleman ||], who showed how to calcu- 
late the tunneling rate by applying the instanton method to a scalar field having 
two nondegenerate minima. According to this method, the nucleation rate is cru- 
cially determined by the value of the action of a nontrivial classical Euclidean 
solution having 0(4) symmetry, called the bounce solution. Coleman also argued 
that the formation and evolution of the true vacuum bubble are described by the 
analytic continuation of the bounce solution to Lorentzian time. Then Coleman 
and De Luccia made an attempt to include gravitational effects in the bounce 
and obtained several important new results. However, their explicit results as- 
sume the thin-wall approximation, under which the space outside the bubble wall 
is exactly in the false vacuum. Hence it was unnecessary to know the behavior of 
the bounce solution there. 

Once the thin-wall approximation is relaxed, it becomes indispensable to con- 
sider the global structure of the bounce solution. In particular, when the Eu- 
clidean Einstein equations are solved simultaneously, a finite vacuum energy den- 
sity in the false vacuum inevitably gives the topological structure of S 4 for the 
bounce solution, leading to the absence of the asymptotic infinity in the Euclidean 
space. Subsequent studies by various people have shown that this fact can cause 
very different behavior of the bounce solution compared with the case of a flat 
Euclidean background. For example, in the other extreme limit in which the 
width of the bubble wall would be very large, the solution reduces to a perfectly 
homogeneous S" 4 with the value of the scalar field staying at the top of the poten- 
tial barrier. This is known as the Hawking-Moss solution ||. Another interesting 
case is that of a potential with degenerate minima. A bounce solution with finite 
action exists in this case because of the compactness of the four-volume, contrary 
to the case of a flat Euclidean background. This describes the formation of a 
topological defect (a domain wall in the case of a single real scalar field) M. 

In the thin-wall and Hawking-Moss limits one can evaluate the bounce action 
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analytically. However, for intermediate cases it must be evaluated numerically. 
Thus, considering the important role played by the bounce action, it is desirable 
and useful to have a formula to evaluate the general bounce action ready to our 
hands. 

In this paper, we consider a real scalar field with a general quartic potential, 
and evaluate numerically the bounce action by varying the parameters. Then, us- 
ing the results, we construct fitting formulae which can be applied to wide ranges 
of the parameter space. This is a generalisation of the procedure of Ref. J7j which 
was only valid in the flat space limit. In Section 2 we introduce the formalism 
and derive our basic equations. In Section 3 we derive analytic results for the 
Hawking-Moss and thin-wall limits. In Section 4 we present our numerical results 
and various analytic fitting formulae. 



2 Method 

2.1 General tunneling formalism 

Following the presciption given by Coleman and De Luccia || , the tunneling rate 
in de Sitter space in the semiclassical limit is given by 

Toce- B/n } (1) 

where B is the bounce action. It is defined to be 

B = S E [M^,r)}-S E [ ( p }, (2) 

where Se is the Euclidean action, r = it is the Euclidean time, 0o is the false 
vacuum value of the scalar field, and </>i(x, r) is the instanton || with minimum 
Euclidean action. An instanton is a solution of the Euclidean equations of mo- 
tion which interpolates from one side of the potential barrier to the other. S The 
instanton with minimum Euclidean action will dominate the tunneling rate. It 
corresponds to the lowest saddle point in the Euclidean action, whose single neg- 
ative mode interpolates between the false and true vacuum states. In flat space it 
can be proved that the instanton has 0(4) symmetry || and it is usually assumed 
that this is also the case in de Sitter space ||. The Euclidean action for a scalar 
field minimally coupled to gravity is 



S E = J d 4 x y/g 
Assuming 0(4) symmetry gives 



16nG 2 



(3) 



ds 2 = + a 2 (0 dQ\ <P = <P{Cj , (4) 



1 The Hawking-Moss instanton is the degenerate case of this in that it sits on top of the 
potential barrier. 
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and 



2.2 A quartic potential 



5 E = 2^/^{-^a(a 2 + l) + 



(5) 



In this paper we limit ourselves to the case of a quartic potential with a local 
minimum (the false vacuum) at = 0, a local maximum (the top of the potential 
barrier) at <p = (f> + , and a global minimum (the true vacuum) at <fi = 0_. The 
potential can be written as 



V(<j>) = V + -m 2 2 - p<\> 6 + A0 4 . 



(6) 



Using the following change of variables we can eliminate two of the five parameters 
from the dynamics: 



*P = -^<P, p = ma, ( = m{,, 



1 8nGV 2Xm' 

Hn = —\ ; , V— — 



m 



and 



V 3 

87rGm 4 87r m 4 



ft 



It will also be useful to introduce the notation 

e — 1 — v . 

Then 



2vr 2 m 2 



Sv, — 



where 



/rfC{-^p(p 2 + l-i/oV)+P 3 



The Euclidean equations of motion are 



$ + 3^ = 11', 
P 

p 2 - 1 + H y = | P 2 (l^ 2 - u 



and the bounce action is 

S = 5b[^(C),p(C)]-5e 



0,-^sinFoC 

-"0 



(7) 
(8) 

(9) 
(10) 

(11) 
(12) 

(13) 
(14) 

(15) 



Note that U(ip) has its local maximum and global minimum at 



(16) 



respectively, and U(ip-) = when v = 1 so that < v < 1. 



2.3 Gravitational back reaction 

The gravitational back reaction will be negligible if k -C 1 except for in the thin- 



wall limit, e<l and Hq <C 1 (see Section |3.2| ). The bounce action for a general 
potential in the thin-wall limit has been given in Ref. 
reduces to 



10 1, and in our case it 



B 



m 2 Ti 2 



1 



2 3 



e-i«) + H 2 



(17) 



We see that the gravitational back reaction is negligible for k <C Hq or k < e. 
Thus for k <C min{l, max{if ) £ }} the gravitational back reaction is negligible and 
we will assume this throughout the rest of the paper. The Euclidean equations of 
motion then become 

p=-^-smH (, (18) 
-"0 

and 

i> + 3H cot HqC ip = U' , (19) 



and the bounce action becomes 



B 



m 2 2tt 



%f 5 °dt sin 3 i7 C 
Hq Jo 



-i) 2 + Z7(V) 



(20) 



The dynamics now depends on only two parameters, v and H , and so a numerical 
treatment becomes tractable. See Figure 1 for a description of the parameter space 
and a sketch of our results. 



3 Analytic Results 

3.1 The Hawking-Moss instant on 

The Euclidean equation of motion, Eq. (|l9|), always has the trivial solution ip = 
This is called the Hawking-Moss instanton. In flat space it has infinite 
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bounce action and plays no role, but because Euclidean de Sitter space has a 
finite volume, 8n 2 /3H 4 , it has a finite bounce action 



m 2 87r 2 



B = l?mi u ™ (21) 

in de Sitter space and is the dominant instanton for sufficiently large Hq. For our 
simple case of a quartic potential the condition under which the Hawking-Moss 
instanton dominates can be obtained by the perturbative analysis below ITlJ , but 
for a general potential this will not be valid because the bounce action depends on 
the global properties of the potential and not just on the local properties around 
the top of the potential barrier. 

Let uj 2 = —U"(ip+) = 2 — and \ = ip — Assume \ <C 1. Then, from 
Eq. ©, 

X + 3# cot#oCx + ^ 2 X^O. (22) 

This equation has a nontrivial regular solution only for uj 2 = n(n + 3)Hq where 
n — 0,1,2, ... . For uj = it is x — constant. For to <C H Q , we expect it still to 
be an approximate solution, in which case B is maximised for \ = 0. Note that 
we must maximise B to obtain a solution that interpolates from one side of the 
potential barrier to the other. For uj = 2H the regular solution is 

X = AcosH (, (23) 

and for to ~ 2H we expect it still to be an approximate solution. Then for A 
sufficiently small, A 2 uj 2 — 4ifg, 

^flK'-if" 2 - 4 ^^}' (24) 

Therefore for uj 2 < 4Hq, B is minimised for A = 0, i.e., the Hawking-Moss 
instanton; while for uj 2 > 4:Hq, B is minimised for A — > oo, i.e., the Hawking- 
Moss instanton is no longer the dominant instanton. Thus the Hawking-Moss 
instanton dominates for 

Hi >\U- h> + ) ■ (25) 



2 V 2 

Now we will calculate the first-order corrections to Eq. ([H]) for < uj 2 —AHq 
1. From Eqs. (O) and (P), 



X + 3H cot H Cx + 4H 2 X = - (uj 2 - AH 2 ) X - 3 (1 - 2z/^+) X 2 + 2u X 3 . (26) 
This has the approximate solution [12] 



X = A cos HqC + 1 2 ^+ A 2 cos 2F C + O (A 3 cos 3H () , (27) 
2 — 3tp-^ ^ ' 
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where 



.4 2 



7(2- 3V> H 



6(1-W>+) 
Substituting into Eq. (BO) then gives 



4^)+o((^-4F f) 



(28) 



5 



m 2 87r 2 



/i 2 3# 4 
For v = this gives 

and for z/ = 1 it gives 



m 2 47r 2 



7(2-3^ + ; 
120 (1 - vfa 



(u, 2 -4tf 2 ) 2 + 



(29) 



/x 2 81B, 4 



o 



63 
20 



1-^(1-4B 2 ) + 



(30) 



m 2 7r 2 



/x 2 12# 4 



7 / 9 \2 



(31) 



3.2 The thin-wall approximation 

From Eqs. (0) and (0), 

W) = ^ 2 (l-^) 2 -^ 4 - 



(32) 



In the thin-wall approximation, e <C 1 so that the true and false vacua are nearly 
degenerate, and H <C 1 so that a large bubble can fit into the Euclidean de Sitter 
space. The instanton can then be divided into three regions 



1. Inside the bubble where 



ip — ip- ~ 1 and U ~ — -e 



2. The bubble wall where 



1 



C~Cw, i>-U' and U~^P 2 (\-^) 2 . 



Therefore 

ip ~ — v^ZTJ ~ — ^(1 — ip) ( therefore r 
3. Outside the bubble where 



1 



1 + e <^ w 



V> = o. 



(33) 



(34) 



(35) 



(36) 
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Putting it all into Eq. ( p0[) then gives 
m 2 2vr 2 



B 



a 2 FT 3 
m 2 2tt 2 



ii 2 FT 3 

and maximising B with respect to £w gives 



1 /"Cw /"l 

-e / c?C sin 3 H ( + sin 3 # Cw / dipip{l-ip) 

2 Jo Jo 

\ih (3 ~ cos ^°^ W + ^ cosS + \ sin3 ^ w 



(37) 
(38) 

(39) 



Note that pw = 0-/Hq) sinf/oCw — Vy^ 2 + ^0 ^ 1- Substituting back into 5 
gives 

„ m 2 7r 2 1 . , 

B = -r- : ! o - ■ (40) 



taniJoC" 



w 



Ho 

e 



" 3 J^TT 2 (s + J^ThI 



Now we will calculate the first-order corrections in e and H to Eq. (f40|) . As 
before, divide the instanton into three parts. 



1. Inside the bubble where 



i/j = ip_ ~l + 2e and C/ ~ + 4e) . 



2. The bubble wall where, from Eqs. (0), and 

4> ~U' -3eip ~U' + 3eV>(l - • 

Therefore 



and 



^~-(l-ie)V(l-V)-2eV, 
^ 2 + £/ ~ (1 - e)ij 2 (l - ij) 2 + ^o/> 2 (5 - 6ip) 



(41) 

(42) 
(43) 
(44) 



Therefore ^ 2 /2 + U = 0{e 2 ) when = 0{e) and ^ 2 /2 + U = -e/2 + 0{e 2 ) 
when if) = 1 + 0(e). Therefore from Eq. ( ]55| ) the boundaries of the bubble 
wall can be taken to be at ( = ± Ine. 



3. Outside the bubble where 



V> = 0. 



(45) 



Putting it all into Eq. ( p0| ) then gives 
m 2 2tt 2 
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// 2 ff 3 







1 /"Cw+lne 

-e(l + 4e / d(sm 3 H ( 
I Jo 



- I d( sin d # C 



(46) 



which from Eqs. ([35|) and fl4*3]), and after noting that we can use Eq. (^) because 
dB/d(w = 0, eventually gives 



B 



m 2 7r 2 



/i 2 3 



e 2 + H 2 \ e + Je 2 + H 2 



3e 2 



3e 3 
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~ 2 ./^^ + 2(e 2 + tf ( 



o, 



In the limit, H <C e, this gives 

B 



m 2 7r 2 



^ 2 12e 3 

while in the limit e <C -f/o it gives 



1 + lOe - -4 + 



5 



m 2 7r 2 



(47) 
(48) 

(49) 



so that we need to calculate to a higher order in Hq in this limit to get the 
first-order correction. We will now do this. 
Set e = 0. Then 



t/ = ^ 2 (i-V0 2 , 

and by symmetry £w = tt/2Hq. Therefore, from Eq. 



(50) 



2U-6H / d( cottfoC^ 2 



2U 



l + 0(H 2 



Therefore 



\^ 2 + U 

Therefore from Eqs. ((2^) and (|35| 



-V2U 



1 + O (H%) 



(51) 
(52) 



B 



rn 



2 2n 2 r 1 



fi 2 H$ Jo 
m 2 7T 2 
~J^3H$ 



dip 



l--H 2 An 2 



^(1-V) + 



1 — - [7T 
2 



2 M 



(53) 



4 Numerical Results 

We solved Eq. ([T^) numerically with v and Hq taking values in the parameter 
space shown in Figure 1, and calculated the corresponding values of B, which is 
defined by 

B = ^4b. (54) 
12/i 2 v ; 
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The results, with the thin-wall divergence, Eq. (pEO|), factored out, are plotted in 
Figure 2. We then fitted various analytic formulae to the results, first for the 
special cases of the straight boundaries of the parameter space and then for the 
full two-dimensional case. 



4.1 One-dimensional fitting 

v = : a cubic potential 

We obtained B = 27.62 in the v = Hq = limit with the accuracy to the given 
figures, and, from Eq. (|30|), B = 256/27 in the Hawking- Moss limit (H = 1/2). 
We compared the numerically calculated B with 

5 fit = 27.62- 72.55# 2 , (55) 

where the coefficients were chosen to give an exact fit at both boundaries of the 
parameter region. In this case the relative error (B — B^ t )/B is found to be 
< 0.6%. A more accurate fit, again exact at both boundaries, is given by 

B = 27 '.62 - limHl - 24.48# 6 . (56) 

Now the relative error becomes < 2.6 x 10 -4 . 

e = 0: degenerate vacua 

The numerical results for this case have already been obtained in Ref. || and 
applied to the nucleation of spherical domain-wall bubbles during inflation. From 
Eq. (|53|) we can see that B — > 4/i^Q as H Q — > 0. Also, from Eq. (|3l|) , we get 
B = 64 in the Hawking- Moss limit (H = l/y/8). Thus the simplest fitting exact 
at both boundaries is 



(57) 



with a relative error < 2%. A better fit, again exact at both boundaries, is given 
by 

(58) 



1 - 1.887F 2 - 3.649F 4 



In this case, the relative error is < 6 x 10 4 . 



H = : flat space 

This case has already been examined in M. In our notation, his results are 



B = — (l + as + be 2 ) , (59) 
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with 

a = 10.052 and 6 = 16.612. (60) 

This fitting formula is not insufficient but it has its largest relative error (~ 0.17%) 
at v — 0. However, the choice of parameters 

a =10.07 and 6 = 16.55, (61) 

is better, making the fitting exact at both boundaries, and giving a relative error 
< 0.14%. 



4.2 Two-dimensional fitting 

Now we try to give a fitting formula valid in the whole range of parameter space. 
Looking at the rather complicated form of Eq. (f|7]), we can guess that it is not 
easy to obtain a simple fitting formula. Figures 3(a) and 3(b) are plots of B\ and 
B 2 which are defined by 



B 



4(1 + 10.07e + 16.55e 2 

5 (e + 5f 
4(1 + 10.07e + 16.55e 2 
6(e + 5) 2 (l + We) 



B 1 



11 

1 + Y e + 3 1 



2 3£ 3 

+ 



2S 2 



Bo 



(62) 
(63) 



where 6 = ye 2 + Hq. From these plots, it seems better to try to fit B 2 instead of 

B or Bi, despite its more complicated factor. Since this factor is written in terms 
of e and 5, we try to fit B 2 with a polynomial of these two parameters. To keep 
the relative error less than about 1% , at least ten terms seem to be necessary. 
We give one example of the fitting: 



B 2 = l + 5>y£W 



(64) 



ij 



with nonzero coefficients 



Ooi = —0.1617, a 3 = — 5.507, 

an = —11.34, a\ 2 — 30.17, au ■ 

a 20 = 12.09, a 21 = -33.24, a 23 

a 30 = 7.728 , and a 32 = —19.34 . 



-21.69 
41.29, 



(65) 



Factoring out the complicated thin-wall factor was necessary to remove the di- 
vergence in the thin-wall limit. Therefore, for small u, B is expected to be fitted 
by a simpler formula. We give the following example of a fitting for v < 0.4: 



B = — (l + 10.07e + 16.55e : 

F 6 \ 



1 + (a + bv)El 



1 + {c + dv + ev 2 )v 2 H$ 



(66) 
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with 

a = -2.635, b = -0.378, c = 43.6 , d=-130, and e = 163 . (67) 
The relative error of this fitting is < 1.4%. 

5 Conclusions 

We have thoroughly investigated the tunneling rate for a scalar field with a general 
quartic potential in de Sitter space in the limit of no gravitational back reaction. 
In this case, the number of parameters of the theory essentially reduces to two. 
We have numerically calculated the tunneling rates over this two-dimensional pa- 
rameter space. Based on analytic formulae which are valid for extreme values of 
the parameters, we have constructed approximate formulae from the numerical re- 
sults with relative errors of < 1%. Since these approximate formulae incorporate 
the effect of the background curvature, they have a wider range of applicability 
than previously known formulae. Hence our results will provide a better under- 
standing of quantum tunneling phenomena and the associated physical processes 
in the early Universe. 
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Figure Captions 
Figure 1: Parameter space 



The dotted line (given by Eq. (|25|)) marks the boundary between the regions 
where the Hawking-Moss (HM) instanton (see Section |3.1|) and our numerically 
calculated instantons dominate. The corner labeled TW is where the thin-wall 
approximation (see Section |3.2|) is valid. The boundaries labeled CP, FS, and DV 



correspond to the limits of a cubic potential, flat space, and degenerate vacua, 
respectively (see Section fil~|) . The solid lines are a contour plot of our numerical 
results (and of the analytic results in the Hawking-Moss region). B is defined in 
Eq. 



Figure 2: Numerical results 

The numerical results expressed in terms of B/B^w, where -Btw is the thin-wall 
bounce action given in Eq. (^0|). 



Figure 3 

A comparison of two possible choices of functions to be fitted by analytic formulae. 
Bi and B 2 are defined in Eqs. ( p2] ) and (§3|), respectively. 
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